Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases.

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Mapping the Complex Genetic Landscape of Human Neurons.

When somatic cells acquire complex karyotypes, they are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons have been limited by relatively small sample sizes. Here, we developed an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We applied this approach to 2,125 frontal cortical neurons from a neurotypical human brain. This approach identified 226 CNV neurons, as well as a class of CNV neurons with complex karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we found that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contained fewer, but longer, genes.

Nitrogen addition and fungal symbiosis alter early dune plant succession.

Anthropogenic nitrogen (N) enrichment can have complex effects on plant communities. In low-nutrient, primary successional systems such as sand dunes, N enrichment may alter the trajectory of plant community assembly or the dominance of foundational, ecosystem-engineering plants. Predicting the consequences of N enrichment may be complicated by plant interactions with microbial symbionts because increases in a limiting resource, such as N, could alter the costs and benefits of symbiosis. To evaluate the direct and interactive effects of microbial symbiosis and N addition on plant succession, we established a long-term field experiment in Michigan, USA, manipulating the presence of the symbiotic fungal endophyte Epichloë amarillans in Ammophila breviligulata, a dominant ecosystem-engineering dune grass species. From 2016 to 2020, we implemented N fertilization treatments (control, low, high) in a subset of the long-term experiment. N addition suppressed the accumulation of plant diversity over time mainly by reducing species richness of colonizing plants. However, this suppression occurred only when the endophyte was present in Ammophila. Although Epichloë enhanced Ammophila tiller density over time, N addition did not strongly interact with Epichloë symbiosis to influence vegetative growth of Ammophila. Instead, N addition directly altered plant community composition by increasing the abundance of efficient colonizers, especially C4 grasses. In conclusion, hidden microbial symbionts can alter the consequences of N enrichment on plant primary succession.

Epichloë Increases Root Fungal Endophyte Richness and Alters Root Fungal Endophyte Composition in a Changing World.

Plants harbor a variety of fungal symbionts both above- and belowground, yet little is known about how these fungi interact within hosts, especially in a world where resource availability is changing due to human activities. Systemic vertically transmitted endophytes such as Epichloë spp. may have particularly strong effects on the diversity and composition of later-colonizing symbionts such as root fungal endophytes, especially in primary successional systems. We made use of a long-term field experiment in Great Lakes sand dunes to test whether Epichloë colonization of the dune-building grass, Ammophila breviligulata, could alter fungal root endophyte species richness or community composition in host plants. We also tested whether nitrogen addition intensified the effects of Epichlöe on the root endophyte community. We found that Epichloë increased richness of root endophytes in Ammophila by 17% overall, but only shifted community composition of root endophytes under nitrogen-enriched conditions. These results indicate that Epichlöe acts as a key species within Ammophila, changing richness and composition of the root mycobiome and integrating above- and belowground mycobiome interactions. Further, effects of Epichloë on root endophyte communities were enhanced by N addition, indicating that this fungal species may become even more important in future environments.

Bimodal Expression Patterns, and Not Viral Burst Sizes, Predict the Effects of Vpr on HIV-1 Proviral Populations in Jurkat Cells.

Integration site landscapes, clonal dynamics, and latency reversal with or without vpr were compared in HIV-1-infected Jurkat cell populations, and the properties of individual clones were defined. Clones differed in fractions of long terminal repeat (LTR)-active daughter cells, with some clones containing few to no LTR-active cells, while almost all cells were LTR active for others. Clones varied over 4 orders of magnitude in virus release per active cell. Proviruses in largely LTR-active clones were closer to preexisting enhancers and promoters than low-LTR-active clones. Unsurprisingly, major vpr+ clones contained fewer LTR-active cells than vpr- clones, and predominant vpr+ proviruses were farther from enhancers and promoters than those in vpr- pools. Distances to these marks among intact proviruses previously reported for antiretroviral therapy (ART)-suppressed patients revealed that patient integration sites were more similar to those in the vpr+ pool than to vpr- integrants. Complementing vpr-defective proviruses with vpr led to the rapid loss of highly LTR-active clones, indicating that the effect of Vpr on proviral populations occurred after integration. However, major clones in the complemented pool and its vpr- parent population did not differ in burst sizes. When the latency reactivation agents prostratin and JQ1 were applied separately or in combination, vpr+ and vpr- population-wide trends were similar, with dual-treatment enhancement being due in part to reactivated clones that did not respond to either drug applied separately. However, the expression signatures of individual clones differed between populations. These observations highlight how Vpr, exerting selective pressure on proviral epigenetic variation, can shape integration site landscapes, proviral expression patterns, and reactivation properties. IMPORTANCE A bedrock assumption in HIV-1 population modeling is that all active cells release the same amount of virus. However, the findings here revealed that when HIV-infected cells expand into clones, each clone differs in virus production. Reasoning that this variation in expression patterns constituted a population of clones from which differing subsets would prevail under differing environmental conditions, the cytotoxic HIV-1 protein Vpr was introduced, and population dynamics and expression properties were compared in the presence and absence of Vpr. The results showed that whereas most clones produced fairly continuous levels of virus in the absence of Vpr, its presence selected for a distinct subset of clones with properties reminiscent of persistent populations in patients, suggesting the possibility that the interclonal variation in expression patterns observed in culture may contribute to proviral persistence in vivo.

Dog10K_Boxer_Tasha_1.0: A Long-Read Assembly of the Dog Reference Genome.

The domestic dog has evolved to be an important biomedical model for studies regarding the genetic basis of disease, morphology and behavior. Genetic studies in the dog have relied on a draft reference genome of a purebred female boxer dog named "Tasha" initially published in 2005. Derived from a Sanger whole genome shotgun sequencing approach coupled with limited clone-based sequencing, the initial assembly and subsequent updates have served as the predominant resource for canine genetics for 15 years. While the initial assembly produced a good-quality draft, as with all assemblies produced at the time, it contained gaps, assembly errors and missing sequences, particularly in GC-rich regions, which are found at many promoters and in the first exons of protein-coding genes. Here, we present Dog10K_Boxer_Tasha_1.0, an improved chromosome-level highly contiguous genome assembly of Tasha created with long-read technologies that increases sequence contiguity >100-fold, closes >23,000 gaps of the CanFam3.1 reference assembly and improves gene annotation by identifying >1200 new protein-coding transcripts. The assembly and annotation are available at NCBI under the accession GCF_000002285.5.

Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes.

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

Comprehensive identification of somatic nucleotide variants in human brain tissue.

BACKGROUND: Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. RESULTS: Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. CONCLUSIONS: This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

An arbuscular mycorrhizal fungus alters switchgrass growth, root architecture, and cell wall chemistry across a soil moisture gradient.

The abiotic environment can dictate the relative costs and benefits of plant-arbuscular mycorrhizal fungi (AMF) symbioses. While the effects of varying light or soil nutrient conditions are well studied, outcomes of plant-AMF interactions along soil moisture gradients are not fully understood. It is predicted that mycorrhizal associations may become parasitic in extreme soil moisture conditions. Under low soil moisture stress, costs of maintaining a mycorrhizal symbiont may outweigh benefits for the host plant, whereas under high soil moisture stress, the host plant may not need the symbiont. In a factorial growth chamber study, we investigated the effects of a plant-arbuscular mycorrhizal fungus symbiosis along a soil moisture gradient on growth, cell wall chemistry, and root architecture of a biofuel crop, Panicum virgatum (switchgrass). Regardless of soil moisture conditions, we found an increase in the number of tillers, number of leaves, root biomass, and amount of cellulose and hemicellulose in response to root colonization by the arbuscular mycorrhizal fungus. The fungus also increased aboveground biomass and changed several root architectural traits, but only under low soil moisture conditions, indicating a reduction in benefits of the mycorrhizal association under high soil moisture. Results from this study indicate that an arbuscular mycorrhizal fungus can increase some key measures of plant growth and cell wall chemistry regardless of soil moisture conditions but is most beneficial in low soil moisture conditions.

Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology.

Long Interspersed Element-1 (LINE-1) retrotransposition contributes to inter- and intra-individual genetic variation and occasionally can lead to human genetic disorders. Various strategies have been developed to identify human-specific LINE-1 (L1Hs) insertions from short-read whole genome sequencing (WGS) data; however, they have limitations in detecting insertions in complex repetitive genomic regions. Here, we developed a computational tool (PALMER) and used it to identify 203 non-reference L1Hs insertions in the NA12878 benchmark genome. Using PacBio long-read sequencing data, we identified L1Hs insertions that were absent in previous short-read studies (90/203). Approximately 81% (73/90) of the L1Hs insertions reside within endogenous LINE-1 sequences in the reference assembly and the analysis of unique breakpoint junction sequences revealed 63% (57/90) of these L1Hs insertions could be genotyped in 1000 Genomes Project sequences. Moreover, we observed that amplification biases encountered in single-cell WGS experiments led to a wide variation in L1Hs insertion detection rates between four individual NA12878 cells; under-amplification limited detection to 32% (65/203) of insertions, whereas over-amplification increased false positive calls. In sum, these data indicate that L1Hs insertions are often missed using standard short-read sequencing approaches and long-read sequencing approaches can significantly improve the detection of L1Hs insertions present in individual genomes.

Stable integrant-specific differences in bimodal HIV-1 expression patterns revealed by high-throughput analysis.

HIV-1 gene expression is regulated by host and viral factors that interact with viral motifs and is influenced by proviral integration sites. Here, expression variation among integrants was followed for hundreds of individual proviral clones within polyclonal populations throughout successive rounds of virus and cultured cell replication, with limited findings using CD4+ cells from donor blood consistent with observations in immortalized cells. Tracking clonal behavior by proviral "zip codes" indicated that mutational inactivation during reverse transcription was rare, while clonal expansion and proviral expression states varied widely. By sorting for provirus expression using a GFP reporter in the nef open reading frame, distinct clone-specific variation in on/off proportions were observed that spanned three orders of magnitude. Tracking GFP phenotypes over time revealed that as cells divided, their progeny alternated between HIV transcriptional activity and non-activity. Despite these phenotypic oscillations, the overall GFP+ population within each clone was remarkably stable, with clones maintaining clone-specific equilibrium mixtures of GFP+ and GFP- cells. Integration sites were analyzed for correlations between genomic features and the epigenetic phenomena described here. Integrants inserted in the sense orientation of genes were more frequently found to be GFP negative than those in the antisense orientation, and clones with high GFP+ proportions were more distal to repressive H3K9me3 peaks than low GFP+ clones. Clones with low frequencies of GFP positivity appeared to expand more rapidly than clones for which most cells were GFP+, even though the tested proviruses were Vpr-. Thus, much of the increase in the GFP- population in these polyclonal pools over time reflected differential clonal expansion. Together, these results underscore the temporal and quantitative variability in HIV-1 gene expression among proviral clones that are conferred in the absence of metabolic or cell-type dependent variability, and shed light on cell-intrinsic layers of regulation that affect HIV-1 population dynamics.

Testing for loss of Epichloë and non-epichloid symbionts under altered rainfall regimes.

PREMISE: Microbial symbionts can buffer plant hosts from environmental change. Therefore, understanding how global change factors alter the associations between hosts and their microbial symbionts may improve predictions of future changes in host population dynamics and microbial diversity. Here, we investigated how one global change factor, precipitation, affected the maintenance or loss of symbiotic fungal endophytes in a C3 grass host. Specifically, we examined the distinct responses of Epichloë (vertically transmitted and systemic) and non-epichloid endophytes (typically horizontally transmitted and localized) by considering (1) how precipitation altered associations with Epichloë and non-epichloid endophytic taxa across host ontogeny, and (2) interactive effects of water availability and Epichloë on early seedling life history stages. METHODS: We manipulated the presence of Epichloë amarillans in American beachgrass (Ammophila breviligulata) in a multiyear field experiment that imposed three precipitation regimes (ambient or ±30% rainfall). In laboratory assays, we investigated the interactive effects of water availability and Epichloë on seed viability and germination. RESULTS: Reduced precipitation decreased the incidence of Epichloë in leaves in the final sampling period, but had no effect on associations with non-epichloid taxa. Epichloë reduced the incidence of non-epichloid endophytes, including systemic p-endophytes, in seeds. Laboratory assays suggested that association with Epichloë is likely maintained, in part, due to increased seed viability and germination regardless of water availability. CONCLUSIONS: Our study empirically demonstrates several pathways for plant symbionts to be lost or maintained across host ontogeny and suggests that reductions in precipitation can drive the loss of a plant's microbial symbionts.

Comparison of village dog and wolf genomes highlights the role of the neural crest in dog domestication.

BACKGROUND: Domesticated from gray wolves between 10 and 40 kya in Eurasia, dogs display a vast array of phenotypes that differ from their ancestors, yet mirror other domesticated animal species, a phenomenon known as the domestication syndrome. Here, we use signatures persisting in dog genomes to identify genes and pathways possibly altered by the selective pressures of domestication. RESULTS: Whole-genome SNP analyses of 43 globally distributed village dogs and 10 wolves differentiated signatures resulting from domestication rather than breed formation. We identified 246 candidate domestication regions containing 10.8 Mb of genome sequence and 429 genes. The regions share haplotypes with ancient dogs, suggesting that the detected signals are not the result of recent selection. Gene enrichments highlight numerous genes linked to neural crest and central nervous system development as well as neurological function. Read depth analysis suggests that copy number variation played a minor role in dog domestication. CONCLUSIONS: Our results identify genes that act early in embryogenesis and can confer phenotypes distinguishing domesticated dogs from wolves, such as tameness, smaller jaws, floppy ears, and diminished craniofacial development as the targets of selection during domestication. These differences reflect the phenotypes of the domestication syndrome, which can be explained by alterations in the migration or activity of neural crest cells during development. We propose that initial selection during early dog domestication was for behavior, a trait influenced by genes which act in the neural crest, which secondarily gave rise to the phenotypes of modern dogs.

Leaf endophytic fungus interacts with precipitation to alter belowground microbial communities in primary successional dunes.

Understanding interactions between above- and belowground components of ecosystems is an important next step in community ecology. These interactions may be fundamental to predicting ecological responses to global change because indirect effects occurring through altered species interactions can outweigh or interact with the direct effects of environmental drivers. In a multiyear field experiment (2010-2015), we tested how experimental addition of a mutualistic leaf endophyte (Epichloë amarillans) associated with American beachgrass (Ammophila breviligulata) interacted with an altered precipitation regime (±30%) to affect the belowground microbial community. Epichloë addition increased host root biomass at the plot scale, but reduced the length of extraradical arbuscular mycorrhizal (AM) fungal hyphae in the soil. Under ambient precipitation alone, the addition of Epichloë increased root biomass per aboveground tiller and reduced the diversity of AM fungi in A. breviligulata roots. Furthermore, with Epichloë added, the diversity of root-associated bacteria declined with higher soil moisture, whereas in its absence, bacterial diversity increased with higher soil moisture. Thus, the aboveground fungal mutualist not only altered the abundance and composition of belowground microbial communities but also affected how belowground communities responded to climate, suggesting that aboveground microbes have potential for cascading influences on community dynamics and ecosystem processes that occur belowground.

Modeling Human Population Separation History Using Physically Phased Genomes.

Phased haplotype sequences are a key component in many population genetic analyses since variation in haplotypes reflects the action of recombination, selection, and changes in population size. In humans, haplotypes are typically estimated from unphased sequence or genotyping data using statistical models applied to large reference panels. To assess the importance of correct haplotype phase on population history inference, we performed fosmid pool sequencing and resolved phased haplotypes of five individuals from diverse African populations (including Yoruba, Esan, Gambia, Maasai, and Mende). We physically phased 98% of heterozygous SNPs into haplotype-resolved blocks, obtaining a block N50 of 1 Mbp. We combined these data with additional phased genomes from San, Mbuti, Gujarati, and Centre de'Etude du Polymorphism Humain European populations and analyzed population size and separation history using the pairwise sequentially Markovian coalescent and multiple sequentially Markovian coalescent models. We find that statistically phased haplotypes yield a more recent split-time estimation compared with experimentally phased haplotypes. To better interpret patterns of cross-population coalescence, we implemented an approximate Bayesian computation approach to estimate population split times and migration rates by fitting the distribution of coalescent times inferred between two haplotypes, one from each population, to a standard isolation-with-migration model. We inferred that the separation between hunter-gatherer populations and other populations happened ∼120-140 KYA, with gene flow continuing until 30-40 KYA; separation between west-African and out-of-African populations happened ∼70-80 KYA; while the separation between Maasai and out-of-African populations happened ∼50 KYA.

Y-Chromosome Structural Diversity in the Bonobo and Chimpanzee Lineages.

The male-specific regions of primate Y-chromosomes (MSY) are enriched for multi-copy genes highly expressed in the testis. These genes are located in large repetitive sequences arranged as palindromes, inverted-, and tandem repeats termed amplicons. In humans, these genes have critical roles in male fertility and are essential for the production of sperm. The structure of human and chimpanzee amplicon sequences show remarkable difference relative to the remainder of the genome, a difference that may be the result of intense selective pressure on male fertility. Four subspecies of common chimpanzees have undergone extended periods of isolation and appear to be in the early process of subspeciation. A recent study found amplicons enriched for testis-expressed genes on the primate X-chromosome the target of hard selective sweeps, and male-fertility genes on the Y-chromosome may also be the targets of selection. However, little is understood about Y-chromosome amplicon diversity within and across chimpanzee populations. Here, we analyze nine common chimpanzee (representing three subspecies: Pan troglodytes schweinfurthii, Pan troglodytes ellioti, and Pan troglodytes verus) and two bonobo (Pan paniscus) male whole-genome sequences to assess Y ampliconic copy-number diversity across the Pan genus. We observe that the copy number of Y chromosome amplicons is variable among chimpanzees and bonobos, and identify several lineage-specific patterns, including variable copy number of azoospermia candidates RBMY and DAZ We detect recurrent switchpoints of copy-number change along the ampliconic tracts across chimpanzee populations, which may be the result of localized genome instability or selective forces.

Resolving complex structural genomic rearrangements using a randomized approach.

Complex chromosomal rearrangements are structural genomic alterations involving multiple instances of deletions, duplications, inversions, or translocations that co-occur either on the same chromosome or represent different overlapping events on homologous chromosomes. We present SVelter, an algorithm that identifies regions of the genome suspected to harbor a complex event and then resolves the structure by iteratively rearranging the local genome structure, in a randomized fashion, with each structure scored against characteristics of the observed sequencing data. SVelter is able to accurately reconstruct complex chromosomal rearrangements when compared to well-characterized genomes that have been deeply sequenced with both short and long reads.

An integrated map of structural variation in 2,504 human genomes.

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

Fungal symbiosis and precipitation alter traits and dune building by the ecosystem engineer, Ammophila breviligulata.

Ecosystem engineer species influence their community and ecosystem by creating or altering the physical structure of habitats. The function of ecosystem engineers is variable and can depend on both abiotic and biotic factors. Here we make use of a primary successional system to evaluate the direct and interactive effects of climate change (precipitation) and fungal endophyte symbiosis on population traits and ecosystem function of the ecosystem engineering grass species, Ammophila breviligulata. We manipulated endophyte presence in A. breviligulata in combination with rain-out shelters and rainfall additions in a factorial field experiment established in 2010 on Lake Michigan sand dunes. We monitored plant traits, survival, growth, and sexual reproduction of A. breviligulata from 2010-2013, and quantified ecosystem engineering as the sand accumulation rate. Presence of the endophyte in A. breviligulata increased vegetative growth by up to 19%, and reduced sexual reproduction by up to 46% across all precipitation treatments. Precipitation was a less significant factor than endophyte colonization for A. breviligulata growth. Reduced precipitation increased average leaf number per tiller but had no other effects on plant traits. Changes in A. breviligulata traits corresponded to increases in sand accumulation in plots with the endophyte as well as in plots with reduced precipitation. Sand accumulation is a key ecosystem function in these primary successional habitats, and so microbial symbiosis in this ecosystem engineer could lead to direct effects on the value of these dune habitats for humans.

The genomic landscape of polymorphic human nuclear mitochondrial insertions.

The transfer of mitochondrial genetic material into the nuclear genomes of eukaryotes is a well-established phenomenon that has been previously limited to the study of static reference genomes. The recent advancement of high throughput sequencing has enabled an expanded exploration into the diversity of polymorphic nuclear mitochondrial insertions (NumtS) within human populations. We have developed an approach to discover and genotype novel Numt insertions using whole genome, paired-end sequencing data. We have applied this method to a thousand individuals in 20 populations from the 1000 Genomes Project and other datasets and identified 141 new sites of Numt insertions, extending our current knowledge of existing NumtS by almost 20%. We find that recent Numt insertions are derived from throughout the mitochondrial genome, including the D-loop, and have integration biases that differ in some respects from previous studies on older, fixed NumtS in the reference genome. We determined the complete inserted sequence for a subset of these events and have identified a number of nearly full-length mitochondrial genome insertions into nuclear chromosomes. We further define their age and origin of insertion and present an analysis of their potential impact to ongoing studies of mitochondrial heteroplasmy and disease.